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^ ■ Abstract 

I The shear viscosity of QED plasma at finite temperature and den- 

' sity is calculated by solving Boltzmann equation with variational ap- 

■ proach. The result shows the small chemical potential enhances the 

viscosity in leading-log order by adding a chemical potential quadratic 
term to the viscosity for the pure temperature environment. 
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^ • 1 Introduction 

: 

^ I A novel state of matter, strongly interacting quark-gluon plasma(sQGP) is 

^ I claimed to be found at Relativistic Heavy Ion Collider at Brookhaven Na- 

tional Laboratory IJJ. The measured V2 was found to reach the hydrodynamic 
limit of an almost perfect fluid with very small viscosity at low transverse 
momentum region. It is desirable to explain this near-perfect fluid behavior 
of sQGP from the theoretical points of view[2j. 

In principle, there are two approaches to calculate transport coefficients. 
One is using the Kubo formulae j2j within the thermal field theory, with 
which people evaluated the shear viscosity via resumming an infinite se- 
ries of ladder diagrams [U EJ IH] . The alternative framework is the kinetics 
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theory [7| IH} 1^1 ITU] . Although the transport equations are hard to solve, the re- 
laxation time approximation(RTA) and variational calculus are two popular 
methods to obtain the transport coefficients. In RTA, people use classic ki- 
netic formulae, but involving the relativistic and quantum effects, to estimate 
the shear viscosity [TTl [T21 HH]- Arnold, Moore and YaffejSlini have studied 
the leading-log contribution as well as the full leading order contribution of 
various transport coefficients of the QCD-like theory at high temperature by 
solving the Boltzmann equation with variational approach. The results in 
the two frameworks are coincident in leading-log order except for some factor 
differences. Some publications also demonstrated that the diagrammatic ex- 
pansion of Kubo formula was actually equivalent to the kinetics calculation 
from the linearized Boltzmann equation if all the possible ladder diagrams 
were resummed in scalar field Jl] and in pure gauge theorvfT^ I16j. In addi- 
tion, one should pay attention to the consistency of the Ward identity with 
the ladder resummation 17j in gauge theory. 

However, most works listed above concentrated on the high temperature 
but vanishing chemical potential except Ref. While actually the net 
baryon number in the central &e ball of heavy-ion collision is not zero rigidly 
though small It makes sense to involve this density effect by introducing 
a chemical potential /i, which is much smaller than the temperature, to study 
how it affects the shear viscosity of the plasma. 

In this paper, we shall try to solve the Boltzmann equation by the varia- 
tional method at high temperature with finite density in QED, following the 
scheme in Ref.[8, for high temperature and zero chemical potential. QED is 
a good toy model for the non-Abelian gauge QCD yet simpler in computa- 
tion. We found the shear viscous coefficient is proportional to T^e^/{\n-) 
and modified by a small factor of (1 + 0.13yU^/T^). 

The paper is arranged as following: in the second section, we will review 
the sketch of solving Boltzmann equation by variational method in the ki- 
netics of transport theory and define the shear viscosity in this framework. 
The associated collision processes on the right hand side of Boltzmann equa- 
tion will be calculated in section 3. And in the fourth section, we use the 
variational method to obtain the shear viscosity. Section 5 is conclusion and 
outlook. 

We use the notation P = {po,p) and p = |p|. The momentum denoted 
by a capital letter is the four-component momentum and the lowercase with 
bold face denotes the three-component momentum. 
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2 Boltzmann equation and viscosity 



Considering a system which shghtly deviates from the equihbrium state by a 
small velocity gradient, one can describe it with the one particle distribution 
which is satisfied the Boltzmann equation 

(1) 

where Vp = p = p/p and F is the external force. In the case of shear 
viscosity, the external field is irrelevant and the time derivation on the left 
hand side may be dropped out due to its higher order contribution in spacial 
gradients^. The right hand side of equation([Q) is the collision term which 
takes the usual form of 

C[/](P) = \ [ \Mip,k;p',k')\'i27rr5^'\P + K-P'-K') (2) 

X {/(p)/(k)[l ± /(p)][l ± /(k)] - /(p')/(k')[l ± /(P')][l ± /(k')]} , 

if only 2^2 elastic collisions are involved. Here p, k, p' and k' denote 
the momenta of the incoming and outgoing particles respectively. The mo- 
mentum space integration is a shorthand for J (2^^^, and \A4\'^ is the 
two-body scattering amplitude. The 1 ± / factor is the final state statistical 
weight for boson with the upper sign and for fermion with the down sign. 

Following Ref . [8 , we expand the distribution function in the near-equilibrium 
state and obtain the linearized Boltzmann equation: 

Sij{p) =Cxij{p), (3) 

where C is the linearized collision operator. By defining the inner product 
in function space and variating the trial function Xij{p), one can obtain the 
shear viscous coefficient 

2 

V '^^^Q max (4) 



where 



Xij{p) = Iij{p)x{p) = \l\{PiPj - \^ij)x{p). (6) 
S,,) = -f^'Y. [ P fo iP)[^ ± /o (P)]X'^(P) (7) 
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Figure 1: The possible processes which contribute to the leading- log in the 
collision term in QED plasma. The solid line is for electron and the wiggly 
line is for photon. 



and the collision term at the right hand side of Boltzmann equation is 

abed 



[ y^\M''^\^ {2TiY 6^*\P + K - P' - K') 
ip,k,k,k' . , 



X 



X 



/o"(p)/o^(k)[l±/oV)][l±/o(k')] 



(8) 



where a, 6, c and d are for species of particles. 

In the above definitions, we adopted the formalisms developed by Arnold, 
Moore and Yaffe|HI with the only differences in the distribution functions 
which involved the chemical potential in the initial and the final states. 
Another notation one should notice is the sum in front of the matrix ele- 
ment which means all possible collision processes relevant to the leading-log 
contribution are involved and properly treated without double counting or 
multi-counting. 



3 Collision terms 

In QED, all possible reactions can be classified as two categories: processes of 
exchanging a boson(Fig.l(a)) and processes of exchanging a fermion(Fig.l(b) 
and (c)), in which the later includes the pair production and the Compton 
scattering processes. Notice that the s-channel scattering is omitted because 
it is infrared finite thus does not contribute to the leading- log result. 

Before going into the next step of calculation, we should specify some 
important approximations and definitions. 



4 



• In our discussion, we adopt the hard forward scattering approximation, 
namely the momentum transfer q ^ eT which is small for all the time 
since it is sufficient to compute the leading-log viscosity. So we neglect 
all the differences between the distribution functions such as /(p) and 
f{p')- The fermion mass is also omitted in this case, for it is in order of 
eT which is much smaller than the hard scale T. Thus the kinematics 
of the two-body collision gives 

cosOpk = I + {I - cos^ 0){1- cos 0) (9) 

where 6pk is the angle between p and k. 6 is the angle between p and 
q, and the angle between k and q as well, since they are approximately 
equal in the forward scattering. is the angle between the p-q plane 
and p-q plane. 

• Due to the energy-momentum conservation, only three of the four mo- 
menta of incoming and outgoing particles are independent. If we prop- 
erly label the particles coming from the same vertex with the similar 
momentum symbols as shown in the Figure 1(a), for example P and 
P', all the three Mandelstam variables can be defined as s = [P + K)^, 
t = {P- P'f and u = {P- K'f. 

• As to the infrared divergence, the two categories of collisions behave 
different. When the momentum transfer q = |p — p'| goes to zero in 
the forward scattering, one finds the infrared singularity in the fermion- 
exchange process is logarithmical while in the boson exchange process 
it is quadratic. Fortunately that is not so bad for the latter case be- 
cause if carefully considering the [x" + ^ X'^ ~ X'^V one may 
find a small emerges which softens the quadratic divergence into 
a logarithmical one. Since now all the collision integrations are log- 
arithmically divergent, the limit cut-offs play important roles in our 
calculation. For transport coefficients like shear viscosity, these inte- 
grations are dominated by the hard scale T of the system which can 
be chosen as the ultraviolet cut-off. As to the infrared limit, the hard 
thermal loop self-energy scale eT is sufficient [8^ . Even in the finite den- 
sity case, the small chemical potential only modifies the infrared cut-off 
by adding a factor like e/i behind eT, which does not contribute to the 
leading- log order In 1/e since we assume the chemical potential is much 
small than the typical momentum scale T, i.e. < eT <^ T. Therefore 
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we will not carefully treat the dq integration and just adopt T and eT 
as the upper and down limits respectively. 

Now let us continue our calculation. 5^(p + k — p' — k') in the integrand 
of equation© helps to perform the k' integration, yet to the 6 function of 
energy conservation, one may introduce a dummy integration variable cuj?] 

/oo 
dujS{uj + p~p')5{uj-k + k'). (10) 
-oo 

With this trick we can integrate over the angles and the remaining integrals 
are 

cq POO POO /•2iT 



POO pq POO POO pZn 

{Xij^Cxij) = j-r^ / dq du dp dk dcpy^lM 
I4vrj Jo J.g Jo Jo Jo ^^^^ 



aft 1 2 
cd I 



X /o(p)/o(k)[l±/o(p)][l±/o(k)] 



X 



[xUp) + - 4-(p') - XU^')T- (11) 

with p' = p + uj and k' = k — uj. For the sake of convenience, we adopt f{p) as 
the fermion distribution and h{p) for boson function in the equilibrium state 
in the following calculation. 



3.1 Boson-exchange Processes 

Unlike the pure temperature case, the system with finite chemical poten- 
tial requires more careful treatment to distinguish the different species of 
fermions with different distribution functions. For the boson-exchange pro- 
cess, Bhabha scattering e^e~ — >■ e^e~ and M0ller scattering e^e^ e^e^ or 
e^e^ — > e^e^ have been involved. Since the s-channel has been omitted since 
it does not contribute to the leading-log order, the distribution functions in 
the Boltzmann equation for both scatterings are 

Bhabha scattering : 2f{p)f{k) [1 - f{p)] [1 - f{k)] 

+2f{p)f{k)[l-f{p)][l-f{k)], (12) 

M0ller scattering : A f{p)f{k)[l - f{p)][l - f{k)] (13) 

where the extra factor of 4 in the M0ller scattering process is from the sum 
over the initial and final states, and the factor of 2 in the Bhabha scat- 
tering comes from the t-channel and u-channel contributions respectively. 
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fip) = ^e^(p+^^) + is the distribution function for the positron and 
f{p) = [e^*^^"^-* + is the distribution function for the electron. 

In the forward scattering, one can easily check s ~ —u, thus the matrix 
element for t-channel is 

8e4^!±^ ~ IQe^'!^ = 16e^^(l - cos0)l (14) 



q 

where the spins of initial and final states have been summed. As to u- 
channel, the matrix element is identical with that of t-channel as long as the 
momentum symbols are well defined. 

In the case of q being small and the particle species a, b being identical 
to c, d{oT d, c) respectively, one finds 

X-,(P')-X-,(P) = q-Vx?,(p) + --- (15) 



^Iij{p)x^{p)' - ^J^C^^PiPi - QtPj - QjPi 



P 

where x^ip)' = dx^{p)/dp. The square of the above equation one obtains 

[X^(p') - X'ip')? = co'irip)? + ^^^^IX^P)? + 0{q'). (16) 

Here, the electron and positron have the same departure from the equilibrium 
which is denoted by x^- One can prove that the cross terms like [xfj(p') ~ 
Xij(p)] ■ [Xi7(k') — Xij(k)] vanishes when carrying out du and dcp integration 
with the factor (1 — cos(/))^ coming from the matrix element. 

Combining the equations (fTT|) and (fT^ - (fT^ and completing the doj and 
d(f) integration, we obtain the collision term for the boson-exchange process 

{X^,,Cx^,t^ = Tdp rdkp'k'{p'[x^{p)f + 6[r{p)r} (17) 

X {f{p)f{k)[l - /(p)][l - fik)] + f{p)f{k)\l - f{p)][l - f{k)] 

+ mfik)[l - f{p)][l - f{k)] + f{p)f{k)\l - f{p)]\l - f{k)]}. 

where we have replaced k with p in the x-functions and placed an extra factor 
of 2 in front of the remaining integration. 
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Noticing that the k-integration can be done after expanding the distribu- 
tion functions in terms of /i/T 



dkk^f{k)[l-f{k)]=T' Y + + 



22^2 ' ^j'Z' 



dkk^f{k)[l-f{k)]=T' 



--ln4^ + 4 + C^(4 

6 T ^T3' 



(18) 
(19) 



we obtain 



Ini 



x(i+^|^) {p'[r{p)r+^[x'{p)n- 



(20) 



3.2 Pair Production 



The pair production process is described by Fig. 1(b) and its reversed process. 
The typical matrix element for this process is: 



ee 1 2 
77 1 



U t 



2u 



t ^ u t 



8e 



2pk^ 
q2 



cos I 



(21) 



Adding up the distribution functions contributions, the equation (|TT|l is recast 
into 



iXijXx 



16e^ 
(4^ 



dq duj j dp 



dk 



2iv 



d(f){l — cos( 



2pk 



{f{p)f{k)[l + h{p)][l + h{k)] + f{p)f{k)[l + h{p)][l + h{k)] 
+b{p)b{k)[l - f{p)][l - f{k)] + b{p)b{k)[l - f{p)][l - m]} 

X IxUp) + 4(k) - XW) - 4(k')]'- (22) 
Expanding the x-function term and ignoring the momenta difference between 
the incoming and outgoing particles we get 

[x-,(p) + X-,(k)-x7,(p')-xZ,(k')f 



+2k,{p) ■ h,{k)[r{p) - X'{P)] [X%k) - x\k)]. (23) 
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And noticing ifAp) = 1 and 



Iij(p) ■ Iij{k) = -{3cos^6pk - 1) = -P2 (cos 6'pfc) 



(24) 



where P2 (cos is the second Legendre polynomial, one can check the cross 
term vanishes when integrating over d(f). We carry out the k-integration by 
expanding the integrand in terms of small fi/T and find, 



dk k f{k)[l + b{k)] = 

POO 

/ dk k f{k)[l + b{k)]=T'^ 
Jo 

POO 

/ dk kbik)[l- f{k)]=T^ 
Jo 

POO 

/ dk k b{k)[l - f{k)] = T"^ 
Jo 



!!:!_ 0.963^ + 0.2984 + (^(4) 



— + 0.963^ + 0.2984 + 0(4' 



0.270^ -0.0484 + (^(4' 

g 2^2 ^J^3' 



— + 0.270^ -0.0484 + (^(4) 
g 2^ 2^2 ^2^3'' 



(25) 
(26) 
(27) 
(28) 



And then the equation ()22|l becomes 
pe^ In -1 



2%5 
+/(p)[l + 6(p)] 



dpp < fip)[l + b{p)] 



TT 



0.963^ + 0.298^^ 



2 2 

— + 0.963- + 0.298^ 
8 T T2 



+6(p)[l-/(p)](y -0.27^ - 0.048^ 



+Kp)[l-/(p)](y + 0.27|- 0.048|^ 



(29) 
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3.3 Compton Scattering 

The Compton scattering process involves both electron and position contri- 
butions. The matrix element for this process is 

\MT/ = -8e'^ = 8e^^(l - cos0). (30) 

The distribution functions for this process is 

f{p)b{k)[i - f{k)][i + b{p)] + f{p)b{k)[i - m][i + b{p)]. (31) 

The x-function terms becomes 

[X^,(p)+X7,(k)-X^,(k)-Xj,(p)]' 

— [rip) - X'{P)T + [r{k) - X^{k)Y. (32) 
After finishing the integration over dk we can recast the equation ljllll into 

/3„4 i„ 1_ i>oo r / 2 2 \ 

ix^^Cx^,P = ^l^ I dpP |/(p)[l + Kp)](y- 0.270|- 0.048|jj 

+ m[l + b{p)] (^^ - 0.270| - 0.048^) I [xUp) ' xlM'- (33) 

4 Variational Method 

As far as shear viscosity is concerned, two species of particles are involved 
and x{p) must take two components 




Accordingly the collision operator C is a 2 x 2 matrix. The left hand side of 
Boltzmann equation© reads 

iX^,,S.,) = -^y^ dpp' [b{p)[l+bip)]x'{p) 

+f{p)[l - f{p)]x\p) + Kp)[1 - Kp)]x\p)]- (35) 
10 



Since we have already obtained all the collision terms in Boltzmann equa- 
tion, we are going to solve the equation 

{Xij,Sij) = {Xij,Cxij) (36) 

to get the shear viscosity by variating the ansatz Xij to reach its maximum 
value. We are not going to argue much about the accuracy of this method in 
this paper, because Arnold et a/jHI have compared it with the exact results at 
high temperature but zero chemical potential environment. And we will see 
the ansatz is the function only in terms of the momentum and the thermal 
variables, thereby we can safely use the same ansatz form in the small /x. 

Before we choose the exact ansatz of and x^? we prefer to demonstrate 
the scheme of this variational calculus. For simplicity all the subscripts and 
momentum dependence of each function and operator are dropped out, and 
the Boltzmann equation becomes 

(x,5) = (x,Cx). (37) 
Expanding the x-function in a finite basis set 

AT 

X = ^ am(pm = a-(j) (38) 

m,=l 

one finds the equation (j?7jl becomes 

am(0m, S) aman{(t)m, C(j)m)- (39) 

m mn 

Redefining S and C in the (pm basis set, one finds 

am C an (40) 



mn 



with S = {(t>m, S) and C = (</>„,, C0„). It's a trivial exercise to give a = C 
and 

V = ^Qma. = YE'^-s = YE^^'^^^'s. (41) 

For the real two-component x-function one can expand it in the same 
finite basis set 

N N 

X'ip) = am(j)m{p), X'(P) = a7V+m0m(p), (42) 

m=l m=l 



11 



where {oi, 02, ■ ■ ■ , a27v} are the independent variational parameters, and adopt 
the one function of the set with natural ansatz (j){p) = p^, 

By using this form of ansatz and neglecting the higher order than /i^/T^, one 
can evaluate 

/?2 3 / b{p)[l+b{p)] 2 



(44) 



S = o dpp'^ \ _ _ ]p 

120^(5)T^ 



\ ^(1 + 0.869^ 



by expanding the fermion distribution function in term of /i/T and neglecting 
the higher order of /i^/T^. 

The collision term C can be obtained likewise by combining equations 

and dnnD 



^ vrT^e^lne 



9 



1(1 + 0.738^) I 128 \-l 1 



(45) 



Inserting equation (j44p and (j45p into equation (j4H) we obtain the shear 
viscous coefficient for QED plasma 

^QED = 187.13-^^ (1 + 0.13^ + 0(^)) , (46) 

which recovers the result of Ref . at /i = and has similar structure as that 
from relaxation time approximation [TH| IT^ 



5 Discussion and Outlook 

So far we have obtained the shear viscosity of QED plasma at finite temper- 
ature and density in the leading-log order. The chemical potential modifies 
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the result in pure temperature case by a small factor of /x^/T^, which ensures 
the modification factor is irrelevant to the sign of net charge of plasma due 
to the symmetry. In addition, the sign in front of the modification factor 
is positive, which indicates that the chemical potential increases the shear 
viscosity of the plasma. Although we obtain this result in small /j, limit, the 
tendency keeps unchanged in the whole region of yU < T. 

In the thermal field theory, we can also obtain such kind of result like 
equation(|mj) by replacing the damping rate by the transport damping ratepUj 
in the boson-exchange case. The reason for this replacement is clear when 
one looks into the kinetic theory: the extra coming from the x-function 
in ee ee scattering softens the quadratic divergence into a logarithmical 
one. This extra small is appeared only in boson-exchange process and is 
the origin of extra sin^ | in the transport damping rate. Carrington, Defu 
and Kobes also pointed out|5j, these x-terms can be explained as an infinite 
series of resummed ladder diagrams. These facts imply that the one-loop 
calculation with usual interaction rate is not complete even in the amplitude 
of order. But the replacement of transport rate improves the calculation and 
makes the results reliable. 

We have calculated the viscosity of plasma involving only 2—^2 processes 
to leading logarithm. But the inelastic scatterings and interference effects 
might be important if we go beyond the leading-log and obtain the complete 
leading order contribution. Furthermore, to explain the near-perfect property 
of QGP, one need to treat the strong coupling system. In this case we have to 
use Kubo formula and calculate the correlation functions of relevant currents. 
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